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We investigate the entanglement transfer from a bipartite continuous-variable (CV) system to a 
pair of localized qubits assuming that each CV mode couples to one qubit via the off-resonance 
Jaynes-Cummings interaction with different interaction times for the two subsystems. First, we 
consider the case of the CV system prepared in a Bell-like superposition and investigate the condi- 
tions for maximum entanglement transfer. Then we analyze the general case of two-mode CV states 
that can be represented by a Schmidt decomposition in the Fock number basis. This class includes 
both Gaussian and non Gaussian CV states, as for example twin-beam (TWB) and pair-coherent 
(TMC, also known as two-mode-coherent) states respectively. Under resonance conditions, equal 
interaction times for both qubits and different initial preparations, we find that the entanglement 
transfer is more efficient for TMC than for TWB states. In the perspective of applications such 
as in cavity QED or with superconducting qubits, we analyze in details the effects of off-resonance 
interactions (detuning) and different interaction times for the two qubits, and discuss conditions to 
preserve the entanglement transfer. 

PACS numbers: 03.67.Mn, 42.50.Pq 



I. INTRODUCTION 



Entanglement is the main resource of quantum infor- 
mation processing (QIP). Indeed, much attention has 
been devoted to generation and manipulation of entan- 
glement either in discrete or in continuous variable (CV) 
systems. Crucial and rewarding steps in the development 
of QIP are now the storage of entanglement in quan- 
tum memories [l|, Q and the transfer of entanglement 
from localized to flying registers and viceversa. Indeed, 
effective protocols for the distribution of entanglement 
would allow one to realize quantum cryptography over 
long distances Q, as well as distributed quantum com- 
putation Q and distributed network for quantum com- 
munication purposes. 

Few schemes have been suggested either to entangle 
localized qubits, e.g. distant atoms or superconduct- 
ing quantum interference devices, using squeezed radi- 
ation [j| or to transfer entanglement between qubits and 
radiation 1,0,1- As a matter of fact, efficient sources 
of entanglement have been developed for CV systems, 
especially by quantum-optical implementations @. In- 
deed, multiphoton states might be optimal when con- 
sidering long distance communication, where they may 
travel through free space or optical fibers, in view of the 
robustness of their entanglement against losses [l(| • 

The entanglement transfer from free propagating light 
to atomic systems has been achieved experimentally in 
the recent years 0, [ll|. From the theoretical point of 
view, the resonant entanglement transfer between a bi- 
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partite continuous variable systems and a pair of qubits 
has been firstly analyzed in [l2| where the CV field is as- 
sumed to be a two-mode squeezed vacuum or twin-beam 
(TWB) state [l3| with the two modes injected into spa- 
tially separate cavities. Two identical atoms, both in the 
ground state, are then assumed to interact resonantly, 
one for each cavity, with the cavity mode field for an 
interaction time shorter than the cavity lifetime. More 
recently, a general approach has been developed [||, in 
which two static qubits are isolated by the real world 
by their own single mode bosonic local environment that 
also rules the interaction of each qubit with an exter- 
nal driving field assumed to be a general broadband two 
mode field. This model may be applied to describe a 
cavity QED setup with two atomic qubits trapped into 
remote cavities. In Ref. [7| the problems related to differ- 
ent interaction times for the two qubits are pointed out, 
either for atomic qubits or in the case of superconducting 
quantum interference devices (SQUID) qubits. The pos- 
sibility to transfer the entanglement of a TWB radiation 
field to SQUIDS has been also investigated in Q. 

Very recently, in |14j the entanglement transfer pro- 
cess between CV and qubit bipartite systems was inves- 
tigated. Their scheme is composed by two atoms placed 
into two spatially separated identical cavities where the 
two modes are injected. They consider resonant inter- 
action of two- mode fields, such as two-photon superpo- 
sitions, entangled coherent states and TWB, discussing 
conditions for maximum entanglement transfer. 

The inverse problem of entanglement reciprocation 
from qubits to continuous variables has been discussed 
in [l5| by means of a model involving two atoms pre- 
pared in a maximally entangled state and then injected 
into two spatially separated cavities each one prepared in 
a coherent state. It was shown that when the atoms leave 
the cavity their entanglement is transferred to the post 
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selected cavity fields. The generated field entanglement 
can be then transferred back to qubits, i.e to another 
couple of atoms flying through the cavities. In a recent 
paper [l6j the relationship between entanglement, mixed- 
ness and energy of two qubits and two mode Gaussian 
quantum states has been analyzed, whereas a strategy to 
enhance the entanglement transfer between TWB states 
and multiple qubits has been suggested in [13]. 

In this paper we investigate the dynamics of a two- 
mode entangled state of radiation coupled to a pair of 
localized qubits via the off-resonance Jaynes-Cummings 
interaction. We focus our attention on the entangle- 
ment transfer from radiation to atomic qubits, though 
our analysis may be employed also to describe the effec- 
tive interaction of radiation with superconducting qubits. 
In particular, compared to previous analysis, we consider 
in details the effects of off-resonance interactions (detun- 
ing) and different interaction times for the two qubits. As 
a carrier of entanglement we consider the general case of 
two-mode states that can be represented by a Schmidt 
decomposition in the Fock number basis. These include 
Gaussian states of radiation like twin-beams, realized by 
nondegenerate parametric amplifiers by means of spon- 
taneous downconversion in nonlinear crystals, as well as 
non Gaussian states, as for example pair-coherent (TMC, 
also known as two-mode coherent) states [l8| , that can be 
obtained either by degenerate Raman processes 11911 or, 
more realistically, by conditional measurements [201 ] and 
nondegenerate parametric oscillators [U [22j. In fact, we 
find that TMC are more effective in transferring entan- 
glement to qubits than TWB states and this opens novel 
perspectives on the use of non Gaussian states in quan- 
tum information processing. 

The paper is organized as follows: in the next Section 
we introduce the Hamiltonian model we are going to an- 
alyze for entanglement transfer, as well as the different 
kind of two-mode CV states that provide the source of 
entanglement. In Section IIIII we consider resonant en- 
tanglement transfer, which is assessed by evaluating the 
entanglement of formation for the reduced density ma- 
trix of the qubits after a given interaction time. In Sec- 
tions IIVI and [V] we analyze in some details the effects 
of detuning and of different interaction times for the two 
qubits. Section IVll closes the paper with some concluding 
remarks. 



II. THE HAMILTONIAN MODEL 

We address the entanglement transfer from a bipar- 
tite CV field to a pair of localized qubits assuming that 
each CV mode couples to one qubit via the off-resonance 
Jaynes-Cummings interaction (as it happens by inject- 
ing the two modes in two separate cavities). We allow 
for different interaction times for the two subsystems and 
assume [3] that the initial state of the two modes is de- 
scribed by a Schmidt decomposition in the Fock number 
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where \nn) = \n) <g> \n) and the complex coefficients 
Cn{%) = (nn\x) satisfy the normalization condition 
S^Lo l c "( :z; )l 2 = 1- The parameter x is a complex vari- 
able that fully characterizes the state of the field. Notice 
that a scheme for the generation of any two-mode cor- 
related photon number states of the form has been 
recently proposed (20| . The simplest example within the 
class ([1]) is given by the Bell-like two-mode superposition 
(TSS): 
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Eq. ([T]) also describes relevant bipartite states, as for 
example TWB and TMC states. In these cases we can 
rewrite the coefficients as c n (x) = co(x)f n (x), where: 



TWB: co(x) = y/1- \x\ 2 f n (x) = x 
TMC : Cq(x) 
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where Io(y) denotes the 0-th modified Bessel function 
of the first kind. For TWB states the parameter |x| 
is related to the squeezing parameter, and ranges from 
(no squeezing) to 1 (infinite squeezing). For TMC 
states |x| is related to the squared field amplitude and 
can take any positive values. The bipartite states de- 
scribed by {T]) show perfect photon number correlations. 
The joint photon number distribution has indeed the sim- 
ple form P n k(x) = S n k\c n (x)\ 2 . For the TSS states the 
joint photon distribution is given by Poo = | Co | 2 and 
P n = 1 - Poo, whereas for TWB and TMC it can be 
written as P n k{%) — 5nkPm{x)\fn{x)\ 2 ■ As we will see 
in the following the photon distribution plays a funda- 
mental role in understanding the entanglement transfer 
process. 

The average number of photon of the states \x), i.e. 
(N)(x) = (x|a^a + wb\x), a and b being the field mode 
operators, is related to the dimensionless parameter |a;| 
by 
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where I\ (y) denotes the 1-th modified Bessel function of 
the first kind. In FigfT] we show the first four terms of 
the photon distribution for TMC and TWB states re- 
spectively as functions of the mean photon number. 

The states in Eq. |T]) are pure states, and therefore we 
can evaluate their entanglement by the Von Neumann en- 
tropy S vn (x) of the reduced density matrix of each sub- 
system. For the TSS case we simply have: 



5 (! 



-Poo log 2 Poo - (1 - Poo) log 2 (l - Poo) (7) 
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FIG. 1: The first four terms P„ n , n = 0, ...,3 of the joint 
photon distribution of the state \x) as a function of average 
photon number (N). (a): TMC, (b): TWB. 



and of course the maximum value of 1 is obtained for 
-Poo = Pn = !• The corresponding state |$+) = |00> J^ 11> 
is a Bell-like maximally entangled state. For TWB states 
the Von Neumann Entropy can be written as: 



S vn (x) = -log 2 (l-\x\ 2 ) 



2|x| 



1 



■log 2 |a;| (8) 



whereas for TMC states we use the general expression: 
S vn (x) = - log 2 Poo(x) - 

oo 

- Poo(x)J2\M*)\ 2 l°E2\fi(x)\ 2 (9) 

i=l 

It is clear that the Von Neumann entropy diverges in the 
limit \x\ — > 1 (for TWB) and in the limit |x| — > oo (for 
TMC) because the probability Poo(x) vanishes. The VN 
entropy at fixed energy (average number of photons of 
the two modes) is maximized by the TWB expression 
©. For this reason TWB states are also referred to as 
maximal entangled states of bipartite CV systems. 

We consider the interaction of each radiation mode 
with a two level atom flying through the cavity. If the 
interaction time is much shorter than the lifetime of the 
cavity mode and the atomic decay rates, we can neglect 
dissipation in system dynamics. On the other hand, we 
consider the general case of atoms with different interac- 
tion times and coupling constants, prepared in superpo- 
sition states, and off-resonance interaction between each 
atom and the relative cavity mode. All these features can 
be quite important in practical implementations such as 
in cavity QED systems with Rydberg atoms and high-Q 
microwave cavities (2^, as noticed in 0- In the inter- 
action picture, the interaction Hamiltonian Hi is given 
by: 



Hi = -HA A (Ja - hA B b ] b + 

+ hg A [a*S%_ + aS?+] + hg B [tfS B 2 
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where S\ 2 ± and S B 2 ± are the lowering and raising atomic 
operators of the two atoms and A a, A b denote the 
detunings between each mode frequency and the corre- 
sponding atomic transition frequency. The initial state 
of the whole system 

|-0(o)> = |s) ® mo)>a ® !#>))* 



evolves by means of the unitary operator U(t) = 
exp[— ifljr] that can be factorized as the product of 
two off-resonance Jaynes- Cummings evolution operators 
Ua(t) and U b (t) [24j related to each atom- mode subsys- 
tem. For the initial state of both atoms we considered 
a general superposition of their excited (|2)) and ground 
(|1)) states: 



A 2 \2) A 
B 2 \2) B 



■Ai\1)a 
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where \A 2 \ 2 + |Ai| 2 = 1 and \B 2 \ 2 + |Bi| 2 = 1. This 
includes the most natural and widely investigated choice 
of both atoms in the ground states, but will also allow us 
to investigate the effect of different interaction times. 

Due to the linearity of evolution operator U(t) and its 
factorized form, the whole system state |V'( r )) a t a time 
t can be written as 

mr)) = 

OO 

= c (x) f n (x)U A (r)\r(0))A ® U B (r)\r(0)) B 
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where \ip n '(0)) a,b = \iP{0))a,b <8 \n)A,B- In each 
of two atom-field subspaces A and B we expand the 
wave-function on the basis {|2)|fc), |l)|fc + 1)}^L U 
{|0)|1)}. The coefficients c A ,i,*(0) -a (2\ A {k\ |V>"(0)>A 
and c B A,k{0) —b (2\ B (k\\%p n (0)) B of the initial states 
are: 



ca,i,o(0) 

CA,2,fe(0) 
CA4,fc+l(0) 



Ai5 n ,o 05,1,0(0) 

A 2 5 k ,n C_B !2 ,fe(0) 
Al8 k +l,n 03,1^+1(0) 



B 2 S k , n (13) 



The Jaynes-Cummings interaction couples only the coef- 
ficients of each variety K whereas ca,i,o(0), cb^^O) do 
not evolve. Therefore, for each variety in the subspaces 
A and B the evolved coefficients can be obtained by ap- 
plying the off-resonance Jaynes-Cummings 2x2 matrix 
Ujk so that: 

C2,fc(r) = l/ii(*,T)o2,fc(0) + l7ia(fc,T)c ll fc + i(0) 
ci, fc+ i(r) = U 21 (k,T)c 2 , k (0) + U 22 (k,T)c hk+1 (0) 



where 

Uu(k,T) 
U 12 (k,T) 
U 22 (k,T) 



,R k r . A . ,R k r 

cos ) — 1 — sin ) 

K 2 ' R k 1 2 ' 



2ig\fk~ 



1 . ( R kT, TJ n . 
- sm (-^-) = U 2 x{k,T), 



Rk 

,R k r A . ,R k r 
cos( — ) + 1^- sm (— ) 



(14) 



where the generalized Rabi frequencies are R k — 
y/4g 2 (k + 1) + A 2 . To derive the evolved atomic den- 
sity operator p\ 2 we first consider the statistical opera- 
tor of the whole system p(r) = \4'( T ))(\' l P( T )\ an d then 
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we trace out the field variables. The explicit expres- 
sions of the density matrix elements in the standard ba- 
sis {\2} A \2} B ,\2) A \1) B ,\1} A \2) B ,\1) A \1} B } are reported 
in Appendix A. 



In the case of both qubits initially in the ground state 
the expression of Af T simply reduces to P22 — 
|/9i4 1, because P22 = P33, and it is possible to derive the 
following simple formula: 



K (Poo,5t) = sin^r) x 



III. ENTANGLEMENT TRANSFER AT 
RESONANCE 

As a first example we consider exact resonance for both 
atom- field interactions, equal coupling constant g and the 
same interaction time r. For the initial atomic states we 
will discuss the following three cases: both atoms in the 
ground state (|1)a|1)b), both atoms in the excited state 
(\2) a \2)b), and one atom in the excited state and the 
other one in the ground state (11)^12)^)- In all these 
cases the atomic density matrix after the interaction p\ 2 
has the following form: 



pI' 2 



/ pu pu 

P 22 

p 33 

\ p* A ,044 



(15) 



The presence of the qubit entanglement can be revealed 
by the Peres-Horodecki criterion [25| based on the exis- 
tence of negative eigenvalues of the partial transpose of 
Eq. (TIB]) . From the expressions of the eigenvalues 
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we see that only X PT can assume negative values. In 
the case of TSS the expression of can allow us to 
derive in a simple way analytical results for the condi- 
tions of maximum entanglement transfer as function of 
dimensionless interaction time gr, as well as to better un- 
derstand the results in the case of TWB and TMC states. 
In order to quantify the amount of the entanglement and, 
in turn, to assess the entanglement transfer we choose to 
adopt the entanglement of formation e p (2|| . We rewrite 
the atomic density matrix in the magic basis [27| p^ B 
and we evaluate the eigenvalues of the non-hermitian ma- 



trix R = p^ {p MBy : 



A l,2 



P22P33 



A 3,4 



U/P17P47±|P14|) 2 (17) 



In this way we calculate the concurrence [28( C = 
max{0, Ai — A 2 — A3 — A4}, where A, = \/X^ are the 
square roots of the eigenvalues X R selected in the de- 
creasing order, and then evaluate the entanglement of 
formation: 



e F = 



l-Vl-C 2 , l-Vl-C 2 
n lo §2 5 



o lo g2 n 



(18) 



x [(1 - P 00 ) cos 2 ( 5 r) - y/(l - P 00 )P( 



mi 



(19) 



We note that only the vacuum Rabi frequency is in- 
volved, a fact that greatly simplifies the analysis of 
atom-field interaction compared to all the other atomic 
configurations. Let us first consider the Bell- like state 
(■Poo = h) ano - look f° r the gr values maximizing the 
entanglement of the two atoms. The solution of equa- 
tion Xf T (^,gT) — —\ is given by gr = f (2fc + 1) with 
k = 0,1,2,.... The above condition is relevant also to 
explain the entanglement transfer for TWB and TMC 
states, as discussed below. To evaluate the entanglement 
transfer also for not maximally entangled TSS states in 
Eq. ([2|), we calculate the entanglement of formation as 
a function of both the dimensionless interaction time gr 
and the probability Poo- As it is apparent from Fig. [2k 
there are large and well defined regions where tp > 0. In 
particular, the absolute maxima {tp = 1) occur exactly 
at Poo = 0.5 and for gr values in agreement with the 
above series. In addition, if we consider the sections at 
these gr values, we obtain exact coincidence with the Von 
Neumann Entropy function S vn (Poo). Therefore, com- 
plete entanglement transfer from the field to the atoms 
is possible not only for the Bell State, though only for 
the Bell State we may obtain the transferral of 1 ebit. 
In Fig. [2t> we consider the entanglement of formation vs 
gr and mean photon number (N) in the TWB case. We 
first note that the regions of maximum entanglement cor- 
respond to those of TSS states and the maxima occur at 
gr values close to f (2fc +1), as shown in Table fl] We 





max 




Poo 


P11 


1.56 


0.87 


0.64 


0.69 


0.21 


4.61 


1.82 


0.81 


0.52 


0.25 


7.85 


1.07 


0.68 


0.65 


0.23 


11.03 


1.07 


0.68 


0.65 


0.23 



TABLE I: Maxima of the qubit entanglement of formation tF 
for the resonant interaction with TWB and for both qubits 
initially in the ground state (see Fig. [2b). 

can explain this by considering the TWB photon distri- 
bution (see Fig. [TJd). We note that the terms Poo and 
Pu are always greater or equal than the other terms P nn 
(n > 1) and that for (N) < 2 they dominate the pho- 
ton distribution (Poo + Pu — 1). Therefore, the main 
contribution to entanglement transfer is obtained from 
the above two terms as for the TSS state. In order to 
explain the absolute maximum found in the second peak 
at (N) — 1.82, we note that in this case Poo and Pu 
are closer to the value 0.5 of a Bell State. In addition, 
for large (N) and gr values, there are small regions (not 
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FIG. 2: Entanglement of formation €f of the qubit systems 
as a function of the dimensionless time gr and the CV state 
parameter Poo (a) or the average number of photons (N) (b,c) 
for the case of both atoms initially in the ground state, a) 
TSS, b) TWB, c) TMC. 



visible in the figure) where entanglement transfer is pos- 
sible. This is due to the terms P nn (n > 1) in the photon 
distribution. In Fig. [2J; we show the TMC case and we 
note that for (N) < 4 there are four well denned peaks 
where the entanglement is higher than in the TWB case. 
Also in this case the gr values of the maxima (see Table 
HI)) nearly correspond to those of TSS states. As in the 





{-^"} max 


£F,max 


Poo Pn 


1.56 


0.89 


0.84 


0.61 0.34 


4.66 


1.09 


0.90 


0.54 0.39 


7.85 


0.99 


0.87 


0.57 0.37 


11.01 


0.99 


0.88 


0.57 0.37 



TABLE II: Maxima of the qubit entanglement of formation 
€f for the resonant interaction with TMC and for both qubits 
initially in the ground state (see Fig. (2ji) . 

previous case this can be explained by the TMC photon 
distribution (see Fig. [T^), where for (N) < 4 the domi- 
nant components of the photon distribution are Poo and 



Pn. The absolute maximum is in the second peak at 
(N) = 1.09 because P 00 and Pn are even closer to the 
Bell State than in the other peaks, and this also explains 
the larger entanglement value cf- For (N) > 4 and large 
gr there are regions with considerable entanglement val- 
ues, due to the fact that Poo and Pn are always smaller 
than the other terms P nn (n > 1) that dominate the 
atom-field interaction. We note that the maxima of ep 
are higher than in the TWB case. 
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FIG. 3: (a): the function A£ T (±,£rr) + ± as in Eq. ((20j) for 
the Bell state case with both atoms in the excited state, (b): 
entanglement of formation cf vs Pqo for TSS states compared 
to Von Neumann entropy (dashed line) for some values of gr, 
corresponding to the numbered minima: (1) 2.03, (2) 4.53, 
(3) 11.07, (4) 26.68. 

A similar analysis can be done in the case of initially 
excited atoms |2)^|2)b- For the TSS states we can again 
write a simple equation for the eigenvalue of the partial 
transpose: 



Af T (Poo, gr) = (1 - Poo) sm 2 (V2 9 T) cos 2 (V2g T ) + 
+ sin 2 ( 5 t)[P 00 cos 2 ( 3 t) - V(l-^oo)Poo cos 2 (\/2 5 t)] 

(20) 



where, with respect to Eq. (|19p . an additional frequency 
is present. For the Bell State we can look for gr values 
maximizing the entanglement transfer. In this case the 
problem can be solved numerically and we found, for ex- 
ample in the range gr = 0—50, that only for gr = 26.68 ~ 
we can solve the equation Af T '(|, gr) — — | with a 
good approximation as shown in Fig. [3^i. In Fig. [3}d we 
consider also non maximally entangled TSS states, show- 
ing the entanglement of formation ep vs the probability 
Poo for gr values corresponding to numbered minima in 
Fig. [SK We see that only for gr = 26.68 a Bell State 
can transfer 1 ebit of entanglement, but the entangle- 
ment transfer is complete also for all the other Poo val- 
ues. A nearly complete transfer can be obtained also 
for gr = 11.07 but in the other cases the entanglement 
transfer is only partial even for the Bell State. We note 
that in [l4| it is shown that for gr — 11.07 one finds 
maximum entanglement transfer for both atomic states 
|2)ji|2)b but starting with a different Bell-like 



_ |10)-|01) 
V2 



field state |*_ 

In the TWB case we find large entanglement transfer 
for gr values very close to the ones of minima (2-4) for 
the TSS states in Fig. [3^. Some gr values correspond- 
ing to maxima of tp in the case of both atoms in the 
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ground state are missing, and the best value of ej? in the 
considered range is at gr = 26.65. Also for the TMC 
states for small (N) we have large entanglement transfer 
corresponding to the above gr values, but in addition for 
(N) > 4 and large gr values there are regions with con- 
siderable entanglement. 

Finally, in the case of one atom in the excited state and 
the other one in the ground state (|1)a|2)b) it is not pos- 
sible to write for TSS states a simple equation as Eq. (|20|) 
because in the atomic density matrix Eq. (fT5")) p 2 2 7^ P33 
unlike in the previous cases. However, there are again 
only two frequencies involved as in the TSS case and we 
can do a similar analysis as for both atoms in the ex- 
cited state. Here we only mention that for dimensionless 
interaction times corresponding to common maxima for 
the different atomic states, the maxima of ep are rather 
lower than in the cases (|1)a|1)b) an d (|2)^|2}b), and 
more in general the transfer of entanglement is sensibly 
reduced as a function of (N) . 

IV. THE DETUNING EFFECT 

^From the practical point of view it is important to 
evaluate the effects of the off-resonant interaction be- 
tween the atoms and their respective cavity fields, which 
can be actually prepared in non degenerate optical para- 
metric processes. We assume equal interaction times for 
both atoms and we consider first the case of resonant in- 
teraction for the atom A and off resonant interaction for 
atom B. As a first example we consider the TMC case, 
both atoms in the ground state and the value gr = 4.66, 
corresponding to the maximum entanglement transfer. 
In Fig. 2k we see that up to detuning values on the or- 
der of the inverse interaction time the entanglement is 
preserved by the off- resonant interaction of atom B. In 
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FIG. 4: The off resonance interaction effect in the TMC case, 
both atoms in the ground state and gr — 4.66. a) Aat = 0, 
and A b t = (dashed line), 1, 2, 3, 4, 5. b) A a t = A s r = 
(dashed line), 1, 2, 3, 4, 5. 

Fig. [4Jd we show the more general case of off-resonance 
for both atoms, taking equal detuning values for simplic- 
ity. The effect is greater than in the previous case but it 
is negligible again up to Agr = 1 . 

Fig.[5]shows the analogous behaviour for the TWB states 
for gr = 4.61. We see that near the peak of entanglement 
the TMC states seem more robust to off resonance inter- 
action than the TWB states. 
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FIG. 5: The off resonance interaction effect in the TWB case, 
both atoms in the ground state and gr — 4.61. a) Aat = 0, 
and A b t = (dashed line), 1, 2, 3, 4, 5. b) A a t = A s r = 
(dashed line), 1, 2, 3, 4, 5. 



V. THE EFFECT OF DIFFERENT 
INTERACTION TIMES 

In the previous analysis we considered equal coupling 
constant and interaction time for both atoms. However, 
experimentally we may realize conditions such that the 
parameter gr is different for the two interactions due to 
the limitations in the control of both atomic velocities 
and injection times or in the values of the coupling con- 
stants 0- 

We first consider the effect of different interaction 
times at exact resonance and simultaneous injection of 
both atoms prepared in the ground state. In Fig.fSkjb we 
show the TMC case for gTA = 4.66 and gTA = 11.01, cor- 
responding to two maxima of entanglement as discussed 
in the previous section, and we investigate the effect of 
different dimensionless interaction times for the atom B 
such that gTB < gTA- We see that increasing the dif- 
ference gTA — gTB the entanglement decreases. However 
for (N) < 4, if gTB has a value close to the one corre- 
sponding to a maximum, as for gTA = 1-56 in Fig. [5k 
and gTA = 7.85 in Fig. [6Jd, i.e. if gTA — gTB — ir, the 
entanglement reaches again large values. The effect is 
more important for gTA = 11.01 where the entanglement 
transfer is the same as for equal interaction times. In 
Fig. [5b, d we show an analogous effect for the TWB case. 

We finally consider the possibility that atom B enters 
the cavity just before atom A. We assume that when 
atom A enters its cavity the two mode field can be still 
described by Eq. |[JJ. Due to the interaction with its 
cavity field the atom B will be in a superposition state 
Bi\1)b + B 2 \2)b- In this case the atomic density ma- 
trix p* 2 has only two null elements, P23 = P32, hence 
we evaluate the eigenvalues of the non-hermitian matrix 
numerically. We calculate the amount of entanglement 
transferred to the atoms after the time r of their simul- 
taneous presence into the respective cavities in the case 
of exact resonance, equal coupling constant and velocity, 
assuming atom A prepared in the ground state. In Fig. [7k 
we show the TMC case for gr — 4.66 and different values 
of \Bi\ 2 ranging from (that is for |1)^|2) B ) to 1 (that 
is |1)a|1)s). We note that the behaviour of gradually 
changes from one limit case to the other one for (N) < 4, 
when the photon distribution approaches that of a TSS 
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FIG. 6: The effect of different interaction times for Aa = 
Ab = 0, and both atoms injected simultaneously in the 
ground state. In a) the TMC case for grA = 4.66 and qtb 
values: a) 4.66, b) 4.4, c) 4.2, d) 4.0, e) 3.8, and 1.56 (dashed 
line). In b) the TMC case for grA ~ 11.01 and grs values: a) 
11.01, b) 10.8, c) 10.6, d) 10.4, e) 10.2 and 7.85 (dashed line). 
In c) the TWB case for grA = 4.61 and grs values: a) 4.61, 

b) 4.4, c) 4.2, d) 4.0, e) 3.8, and 1.56 (dashed line). In d) the 
TWB case for gr A = 11.03 and gr B values: a) 11.03, b) 10.8, 

c) 10.6, d) 10.4, e) 10.2 and 7.85 (dashed line). 

state. For larger values of (N) we see the occurrence of a 
second peak in the case |1),i|2)b. In Fig. |7Jd we show the 
TWB case for gr = 4.61 and we note that the gradual 
change described above occurs for nearly all values of the 
mean photon number (N) . 



a) 




VI. CONCLUSIONS 

In this paper we have addressed the transfer of entan- 
glement from a bipartite state of a continuous-variable 
system to a pair of localized qubits. We have assumed 
that each CV mode couples to one qubit via the Jaynes- 
Cummings interaction and have taken into account the 
degrading effects of detuning and of different interaction 
times for the two subsystems. The transfer of entangle- 
ment has been assessed by tracing out the field degrees 
of freedom after the interaction, and then evaluating the 
entanglement of formation of the reduced atomic density 
matrix. 

We found that CV states initially prepared in a two- 
state superposition are the most efficient in transferring 
entanglement to qubits with Bell-like states able to trans- 
fer a full ebit of entanglement. We have then considered 
multiphoton preparation as TWB and TMC states and 
found that there are large and well defined regions of in- 
teraction parameters where the transfer of entanglement 
is effective. At fixed energy (average number of photons) 
TMC states are more effective in transferring entangle- 
ment than TWB states. We have also found that the 
entanglement transfer is robust against the fluctuations 
of interaction times and is not dramatically affected by 
detuning. This kind of robustness is enhanced for the 
transfer of entanglement from non Gaussian states as 
TMC states. 

Overall, we conclude that the scheme analyzed in this 
paper is a reliable and robust mechanism for the engi- 
neering of the entanglement between two atomic qubits 
and that bipartite non Gaussian states are promising re- 
sources in order to optimize this protocol. Finally, we 
mention that our analysis may also be employed to as- 
sess the entanglement transfer from radiation to super- 
conducting qubits. 



FIG. 7: The entanglement of formation ep vs (N) for a time 
gr of simultaneous presence of both atoms after the delayed 
injection of atom A. Atom A is prepared in the ground state 
and atom B in superposition states such that |-Bi| 2 = 
(dashed line), 0.25 (a), 0.50 (b), 0.75 (c), 1 (dotted line), 
a) the TMC case with gr = 4.66. b) the TWB case with 
gr = 4.61. 
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APPENDIX A: ATOMIC DENSITY MATRIX ELEMENTS 

The elements of the 4x4 atomic density matrix p 12 in the standard basis {12)^12)^, |2}^ 1 1) s , 1 1}^ |2) s , 1 1) ^4 1 1) s } 
are listed below. 



pnW = |co(a ; )| 2 {|A 2 | 2 | J B 2 | 2 £|/,( a; )| 2 |C/ A11 (i,T)| 2 |[/Bii(j,r)| 2 + 



3=0 

OO 

+A\A 2 BIB 2 fj (x)f*+i WUauU, r)U B11 (j, r)U* A12 (j, r)U* m2 (j, r) + 
3=0 

OO OO 

+ |^ 2 | 2 | J B 1 | 2 ^ |. J T 7+1 (^)| 2 |^ 11 0- + 1,^)| 2 |^ 12 0% ^)| 2 + l^ipl^l 2 ^ |/ J + 1 (^)| 2 |^ 12 0% ^)| 2 |^i 1 (.7 + 1,^)| 2 + 

3=0 j=0 
00 

+ |A 1 | 2 |i? 1 | 2 ^|.f, + 1 (x)| 2 |^ 12 (j,T)| 2 |C/ B12 ( J ,T)| 2 + 

3=0 

00 

+A^B 1 B^f3+i^)f;^)UAi2(j,T)U B12 (j,r)UX 11 (j,r)U* B11 (j, tau)} (Al) 
3=0 



P22(x) = |c (x)| 2 {|A 2 | 2 |i? 2 | 2 £ Ifi-^PAiiU l,r)| 2 |f/ S21 (j - l,r)| 2 + 



3 = 1 



+A* 1 A 2 B* 1 B 2 J2f j _ 1 (x)f*(x)U A11 (j - l,r)U B21 (j - l,r)U* A12 (j - l,T)U* B22 (j - l,r) + 



+ \A 2 \ 2 \B 1 \ 2 [J2 \Mx)\ 2 \U a11 (3,t)\ 2 \U B22 (j l,r)| 2 + |^n(0, r)| 2 ] + 



3 = 1 



+ |,4 1 | 2 |i? 2 | 2 £ \fj(x)\ 2 \U A12 (j l,r)\ 2 \U B21 (j,r)\ 2 + |^i| 2 |Bi| 2 £ \fj(x)\ 2 \U A1 2(j - 1, r)\ 2 \U B22 (j l,r)| 2 + 



OO 



+A 1 ^B 1 B 2 *^/ j ( a; )/;_ 1 ( a; )f/ A12 (j - l,r)C/ B22 (j - l,r)I^ n (j - l,r)U* B21 (j - l,r)} (A2) 



p 33 W = |co(a ; )| 2 {|A 2 | 2 | J B 2 | 2 £|/,.(a ; )| 2 |[/ A21 (i,T)| 2 |[/Bii(i,r)| 2 + 



J=0 



+^A 2 i? 1 * J B 2 £ fi{x)fi +1 {x)UA2i{j, T)U B11 (j, r)U* A22 (j, r)U* m2 (j, r) + 



+ |A 1 | 2 |i? 2 | 2 [£ l/.^H^O- - l,r)\ 2 \U B11 (j,r)\ 2 + |^n(0,r)| 2 ] + 



3 = 1 

OO 



+ \A 2 \ 2 \B 1 \ 2 J2 \f j+ i(x)\ 2 \U A21 (j + l,r)| 2 |[/ B12 (j,r)| 2 + {A^B^ 2 ^ l/,+i WH^ 22 (j, r)| 2 |^i 2 (j, r)| 2 + 
j=o i=o 



00 



fiMx)fj{x)UA22(j, r)U Bl2 { 3 , r)U* A21 (j, r)U* B11 (j, r)} (A3) 
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p 4 4(a ; ) = |co(a ; )| 2 {|A 2 | 2 | J B 2 | 2 ^|/ j _ 1 (a ; )| 2 |^ 21 (j-l,T)| 2 |C/B2i(j-l,T)| 2 + 

oo 

+A\A 2 B* 1 B 2 Y J fj-r{x)f*{x)U A21 (j - l,r)U B21 (j - l,r)U A22 (j - l,r)U* B22 (j - l,r) + 
i=i 

OO 

+ |A 2 | 2 |i? 1 | 2 [^|/,(x)| 2 |C/ A21 ( J ,r)| 2 |C/ B22 0--l,r)| 2 + |^ 21 (0,r)| 2 ] + 

OO 

+ |A 1 | 2 |i? 2 | 2 [^|/,(x)| 2 |C/^ 22 (. ? -l,r)| 2 |^2i(j,r)| 2 + |C/ B21 (0,T)| 2 ] + 

OO 

+ |A 1 | 2 |B 1 | 2 []T |/,(x)| 2 |C/^ 22 (. ? - l,r)| 2 + |t/ S22 (j - l,r)| 2 + 1] + 

OO 

+A 1 A*s 1 s*^/ j ( a; )/;_ 1 ( a; )f/ A22 (j - i,r)u B22 (j - i,r)u* A21 (j - i,T)i^ 21 (j - i,t)} (A4) 



+ |A 1 | 2 B 2J B* ]T |.f,(x)| 2 ;7 A12 0- - 1,t)| 2 |C/ b11 (j,t)C/^ 22 0- - l,r) + 



+ J 4 1 A*| J B 2 | 2 ^/ J (x)/;, 1 (x)C/ A12 (j - l.rJI/BnO'.rJI^nO" - 1,t)E£ 2 iO' - 1,t) + 



+ A, A* 2 1 B! | 2 (a;)/; (x) U A12 (j, r) C/ B12 (j, r) C^ n (j, r) C/ B22 (j - 1 , r) + A (x);7 A12 (0, t)(7bi 2 (0, r)C/l n (0, r)]} 

(A5) 



3=0 



+|A 1 | 2 B 1 B 2 *[]T/ j+1 (z)/^ 



i=i 



+AtA 2 |B 2 | 2 [g|/ j (*)| 2 ^ n (j^ 



OO 



+A{A 2 \B 1 \ 2 J2 \f j+ i(x)\ 2 U A11 (j + l,r)\U Bl2 {j,T)\ 2 W A22 { 3 ,r)} 

3=0 



(A6) 
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oo 

p 14 (x) = \c (x)\^\A 2 \ 2 \B 2 \ 2 Y,fj^)f;^)UAii(j,r)U B n(j,r)U* A21 (j - hr)U B2l {j - l,r) + 

OO 

+A\A 2 B* l B 2 [Y J l/,(^)| 2 t/Aii(j,r);7Bii(j,T)C/l 22 (. 7 - l,r)U* B22 (j - l,r) + f/ A11 (0, r)tf B n(0, r)] + 
i=i 

OO 

+ | J 4 2 | 2 |B 1 | 2 [^/ 3 (x).f;_ 1 (x);7 A11 (j + l,T)U m2 (j,T)U* A21 (j,r)U B22 (j - l,r) + 
+/i(a;)C^ii(l,r)^ B i2(0,r)Ca 21 (0,r)] + 

OO 

+ | J 4 1 | 2 | J B 2 | 2 [^/ J+1 (x).f;(x)C/ A12 (j,r)C/ B11 0- + l,r)f/l 22 (j - 1, r)f/ S21 (i, r) + 
+/i(x)[/ a12 (0,t)[/ b11 (1,t)[/ b12 (0,t)] + 

OO 

+ | j 4 1 | 2 | j B 1 | 2 [^/ j+1 (x)/;(x)[/ a12 (j,t)C/ b12 (j,t)C/1 22 (j - l,r)f/ S22 (j - l,r) + 
i=i 

+ fl(x)U A 12(0,T)U B12 (0,T)} + 

OO 

+A 1 ^B 1 B 2 *^/ j+1 ( a; )/;_ 1 ( a; )f/ A12 (i,T)f/ sl2 (i,T)f/l 21 (j - l,r)[/ B21 (j - l,r)} (A7) 



OO 

p 23 (x) - IcoWl 2 ^^^^*^ |/,(*)| 2 ^ii(j,t)[/ B22 (. ? - l,r)U* A22 (j - l,T)U* mi (j,r) + U a11 (0,t)U b11 (0,t)} 

3=0 

(A8) 

OO 

Pu{x) = \c Q {x)\ 2 {\A 2 \ 2 B 1 B* 2 Y J f 3 ( x )fU^ x ) U ^^)UB22{] - l,r)U* A21 (j - l,r)t/* 21 (j - l,r) + 

OO 

+ |A 1 | 2 J B lJ B 2 *[^/ J+1 (x)/;(x)C/ A12 0-,r)C/ B22 (j,r);71 22 0- - 1, r)t/ S21 (j, r) + 
+f 1 (x)U A12 0,r)U B22 (0,r)U B21 (0,T)} + 

OO 

+AIA 2 \B 2 \ 2 [Y / \M^\ 2 Uaii(j,t)\U B21 (j,t)\ 2 UX 22 (j 
i=i 

OO 

+AlA 2 \B 1 \ 2 [Y J \f ] {x)\ 2 U All { 3 ,T)\U B22 { 3 -^ (A9) 



OO 

M*) = IcoCx)! 2 !!^! 2 ^*^^ |/ J (^)| 2 |^A2i(i,^)| 2 C/i3ii(j,^)^ 22 C7 - + |^A2i(0,x)| 2 t/ sll (0,x) 

OO 

+ lAip^i?,^ |/,(x)| 2 C/ A22 (. ? - l,r)U ml (j,r)U* A22 (j - l,r)U* B22 (j - l,r) + ^ B n(0,r)] 

OO 

+ A x A* I i? 2 1 2 ^ (^) x (^) t/^ 2 2 (i - l,r)f/ sll (i,r)f/l 21 (j - l,r)U* B21 (j - l,r) 

OO 

(A10) 
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